d-wave Quasiparticles in the Tilted Vortex Lattice 
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We investigate the quasiparticle spectrum of a d-wave superconductor in the mixed state, fo- 
cussing on the effects of varying the in-plane angular alignment (or tilting) of the vortex and crystal 
lattices. Our approach starts with a linearized Bogoliubov-de Gennes equation cast into a simple 
form by the Franz-Tesanovic singular gauge transformation. The zero magnetic field spectrum is 
characterized by the anisotropy an near the four Fermi-surface nodes; at large etc the model simpli- 
fies into a one-dimensional form that provides asymptotic results. Numerical and analytical results 
are presented and compared with previous studies. We find that tilting can be understood in terms 
of the one-dimensional model at large old- We see a striking dependence on the choice of singular 
gauge, consistent with earlier results, but do not attempt to resolve the issue of which gauge best 
describes the microscopic physics. 



The low-energy quasiparticle spectrum of a d-wave su- 
perconductor in the mixed state has been the focus of 
recent investigations. As the cuprate materials are be- 
lieved to have an energy gap with predominately d x 2_ y 2 
symmetry 1,2 , such work may prove important for under- 
standing certain probes of low-temperature thermody- 
namics and transport in a magnetic field. The recent 
studies have generally restricted their attention to one 
or two fixed vortex lattice geometries, with or without 
small distortions to break the Bravais lattice symmetry. 
However, scanning tunneling microscopy experiments 3 
on YBa2Cu307_5 (YBCO) in a magnetic field of 6 Tesla 
at 4.2 Kelvin show an oblique lattice with an angle be- 
tween primitive vectors of about 77° . This suggests that 
it may be important to consider a range of vortex lat- 
tice geometries if a connection between experiment and 
the calculated spectra is eventually to be made. In this 
paper we investigate the low-energy spectrum as the in- 
plane angular alignment, or tilting, of a square vortex 
lattice and the underlying crystal is varied (Fig. 1). 

In zero magnetic field there are four gapless points on 
the Fermi surface of a single Cu-oxide plane, and the 
spectrum linearized in momentum about each of these 
points takes the form of an anisotropic Dirac cone. Using 
(x',y') to denote crystal lattice coordinates, throughout 
this paper we specialize to a d x 'y' order parameter and 
the node at (0, /c_f), about which the linearized spectrum 
is: 



E(k) = ± 



(1) 



Here va = Ao/pf, where Ao is the maximum magnitude 
of the superconducting gap on the Fermi surface, and 
H = 1. The anisotropy of the spectrum is characterized 
by the ratio ud = vf/va- Experimentally old is known 
to be about 14 for YBCO and 20 for Bi 2 Sr 2 CaCu 2 8 4 . 
We study the mixed state spectrum within the framework 
of the Bogoliubov-de Gennes (BdG) equation 5 , general- 



ized to incorporate a d-wave order parameter and lin- 
earized in momentum operators about the gap nodes 6-9 . 
The Franz-Tesanovic (FT) singular gauge transformation 
can then be used to simplify the problem by casting it in 
the form of an anisotropic two-dimensional (2d) mass- 
less Dirac equation with a periodic vector and scalar 
potential 10 . We will be particularly interested in the 
spectrum at large values of a d , where a one-dimensional 
(Id) model provides asymptotic analytical results 11-13 
that can be compared with the full 2d model. 

We extend the 2D model of Franz and Tesanovic 10 to 
treat tilting at arbitrary angles. Further details of this 
model can be found in the studies by Marinelli, Halpcrin 
and Simon 14 and Vafek et. al. 9 ; here we only sum- 
marize the main features. The magnetic field satisfies 
H c i <C H <C H C 2, implying ( < d < A, where £ is the 
Ginzburg-Landau (GL) coherence length, d is on the or- 
der of the nearest-neighbor vortex separation and A is 
the magnetic penetration depth. We remove the length 
scales A and £ from the problem by taking the magnetic 
field and the magnitude of the order parameter to be spa- 
tially uniform, determining the order parameter phase by 
a linear GL equation. We do not feed the resulting quasi- 
particle spectrum back into the BdG equation to find a 
self-consistent solution, which is a good approximation 
at low temperatures. 

The heart of the FT model is a singular gauge trans- 
formation that requires the decomposition of the vortex 
lattice into A and B sublattices, so that each unit cell 
contains equal numbers of A and B vortices. The spec- 
trum should be independent of the singular gauge used, 
but several studies indicate this is not the case. This was 
first pointed out by Vafek et. al. 9 and discussed further 
in several subsequent works 13,15 ' 16 . We believe this issue 
may be resolved by implementing appropriate boundary 
conditions (b.c.) at the vortex cores, where the FT gauge 
transformation is not valid due to the sharp spatial vari- 
ation of the order parameter magnitude. 
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FIG. 1. Schematic view of the intersection of the flux line 
(vortex) lattice with a single Cu-oxide plane. The magnetic 
field is perpendicular to the plane, which has an (approx- 
imately) square crystal lattice with its orientation described 
by the x' and y' axes. The crystal axes are tilted by an angle 6 
from the square lattice of vortex cores, which are represented 
by the solid circles. Previous work has generally restricted 
attention to 9 = 0° or 45°. 

An analysis of the physics inside the cores will likely be 
necessary to fix the b.c. appropriately. 

We do not attempt to address this issue here, but 
rather study the unregularized theory in two different 
singular gauges, each corresponding to a different choice 
of sublattices. One will be called ABAB (Fig. 2a), the 
other AABB (Fig. 2b). We work in units where d, 
the side of the ABAB unit cell, is unity. This results 
in the following linearized Hamiltonian before singular 
gauge transformation: 



T^BdG 



vf(p v > 

p F \Px' 



{Px',A(r)} 
-v F (p y > + ~ c Ay') 



PF 



(2) 



Here A(r) gives the spatial variation of the GL order 
parameter in the magnetic unit cell unit cell, which is 
aligned with the x and y axes (Fig. 2). Upon making 
the FT transformation and exploiting Bloch's theorem 
to write an effective Hamiltonian at every point (k x ,k y ) 
of the first vortex lattice Brillouin zone, one obtains 



■Ay>)a 3 



(3) 



H iPx' + k x 

a D 



+ A x >)a 1 



A = %(vf 



vf) is the 

B 



where we have set vf = 1. 
effective vector potential and $ — , , , , 
fective scalar potential. vf' B is the superfluid velocity 
associated with the A, B sublattice. 



v~ yl ) the ef- 




FIG. 2. The magnetic unit cells for two choices of singu- 
lar gauge. The square ABAB unit cell (a) has side d and 
is rotated from the underlying vortex lattice by 45°. The 
rectangular AABB unit cell (b) has long side \/2d and is 
aligned with the vortex lattice. In both cases the distance 
between nearest-neighbor vortices is d/\/2, and the tilting of 
the (x',y') crystal coordinates is measured from the vortex 
lattice as shown. The (x,y) coordinate system is aligned with 
the magnetic unit cell. 

Many of the results of Rcf. 14 still hold for this tilted 
Hamiltonian. In particular, the zero energy degenerate 
doublet at the Brillouin zone center (r point) is pre- 
served to all orders in perturbation theory, since this re- 
sult requires only particle-hole symmetry and a Bravais 
vortex lattice 14 . Also, except for special tilt angles and 
anisotropics, we expect no zeroes at other points in the 
zone. This is because tilting spoils the high symmetry 
status of the k y = axis where it was found there could 
be zeroes 14 . However, invariance under k — > — k is pre- 
served. We can therefore restrict our attention to the T 
point to understand the very-low-energy band structure 
and density of states (DOS). 

As pointed out by Vishwanath 15 , the splitting in the 
band structure near any degenerate doublet can be shown 
by first-order perturbation theory to be an anisotropic 
Dirac cone to linear order in the fc-space displacement. 
This actually gives the exact result for the slope of the 
lowest-energy bands as the T point is approached from a 
particular direction. Writing k = fc(cos <j>, sin <f>), we have 



H eff (k) = H eff (0) + fc[sin(0 - 6)<r 3 + — cos(0 - <9)<7 3 ]. 

an 

(4) 

An elementary degenerate perturbation theory calcula- 
tion gives the perturbed energies near the T point: 



e± = ±k\Jh sm 2 ((f> -6)+t 2 cos 2 ((f> - 6) + t 3 sin (2(0 - 6)) + 0(k 2 ), (5) 



where the coefficients can be expressed in terms of ma- 
trix elements of g\ and a 3 between the states of the de- 
generate doublet. Except for the orientation of the new 
anisotropy axes, the anisotropic Dirac cone described by 
Eq. 5 is completely characterized by the maximum and 



minimum values ofe + /k, which we denote by v + and 
respectively. We will study these velocities as a function 
of cxd and 9. 

In order to apply the Id model to obtain informa- 
tion at large an, we must specialize to rational tilt an- 
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gles 9 = arctan(p/q) for integer p and q. Details of 
this model are presented in the papers of Knapp, Kallin 
and Bcrlinsky 12 , and Marinelli and Halperin 13 . The es- 
sential idea is that as ao — > oo the quasiparticles see 
vortices that are effectively smeared out in the "hard" 
^-direction, so that all gauge-invariant quantities be- 
come functions only of the coordinate in the UA-direction. 
Choosing coordinates x and y for the soft and hard di- 
rections, respectively, the Id Hamiltonian for crystal mo- 
mentum k is: 



Hu(k) 



^(p x + k x ) -k y + U B (x)J- 



(6) 



Here U A ' B is the potential encapsulating the residual ef- 
fects of the A, B vortices. For the gauges considered we 
can always choose a unit cell with only two columns of 
smeared-out vortices, provided 9 is rational. This is true 
because the perpendicular distance from one column of 
A/B vortices to the nearest A/B column is always the 
same. 

There are essentially two distinct possibilities when we 
go to the Id model, which depend on the gauge choice and 
on 9. These are illustrated in Figure 3. Configurations 
with distinct columns of A and B vortices (Fig. 3a) were 
shown 13 to give rise to band structure features exponen- 
tially small in old- In particular, for the ABAB gauge 
with 9 = 45°, it was found that t>_ ~ exp(— ao^ 2 
This can easily be generalized when the distance between 
two adjacent columns (of different types) is always the 
same, by scaling the Id unit cell by a factor of l/d v 
in the x-direction and d v in the y-direction, where d v 
is the distance between adjacent vortices in the same 
column. Following the calculation of Ref. 13 one finds 
w_ ~ exp(— a£>TT 2 /Wdf,). 

When the A and B columns overlap (Fig. 3b) the 
asymptotic behavior is very different. It has been ar- 
gued and verified numerically that there are in this case 
no exponentially small features in the band structure 13 . 
We can go further to calculate the exact ID band struc- 



ture on the k v 

j A _ tjB _ 



line, and hence also V-. Since 



U = U = U (x) in this case, we have 

H ld (k x ,0) = U{x) + —{ Px + k x )a 1 . (7) 
a D 

The matrix structure of this Hamiltonian is trivial, and it 
separates into two uncoupled components. The resulting 
differential equations are: 

U{x)4,±— {-i-^ + k x )i> = E±i>. (8) 

Upon integration one finds 

ip{x) = ip e ±iaD r{E±-u(x) T k lc / aD )_ ( Q ) 

ip is a Bloch function, so we must impose periodic bound- 
ary conditions, thus obtaining a quantization condition 



(a) 



(b) 



FIG. 3. Diagram of the two qualitatively different configu- 
rations of the Id model. Dark circles/lines represent one type 
of vortex/column of vortices. The shaded circles and lines 
represent the other. In (a) the A and B vortices form distinct 
columns, while in (b) all columns contain equal numbers of A 
and B vortices. 



for the energies: 

rl/2d v 



/ 1 ' 2d " k 9-rrn 

dx{E ± -U{x)T—) = —. 
-l/2d v a D OL D 



(10) 



For those solutions with \E\ greater than the maximum 
of U(x), this is exactly the Bohr-Sommcrfeld quantiza- 
tion condition used in Ref. 13 to obtain the spectrum 
in a WKB approximation 13 , only here it is exact. At 
the r point (k x = 0), there are doubly degenerate so- 
lutions E± = E^. The degeneracy is split for finite k x 
and E± = E„ ± k x /ao (see Figs. 10 and 11 in Ref. 
13). The lowest bands are simply E± = ±k x /ao, im- 
plying i>_ = l/a_D in striking contrast to the case of 
separated columns. Still more striking is the fact that 
we can switch between exponential (separated columns) 
and algebraic (overlapping columns) decrease of t>_ sim- 
ply by choosing a different singular gauge. For example, 
if the gauge choice did not matter the AABB lattice with 
9 = 45° would be the same as ABAB and 9 = 0°. But 
the spectra are not the same, since the former case has 
overlapping columns in the ID model and the latter does 
not. These conclusions are in excellent agreement with 
the numerical results 13 . 

We have numerically calculated v± as a function of 
an and 9 in the 2d model with a reciprocal space cutoff. 
Working in momentum space allows tilting to be handled 
very simply, and avoids the complications associated with 
discrctizing a Dirac equation on a real space lattice 14,17 . 
The major disadvantage is that the matrices to be diag- 
onalized are not sparse. For the square unit cell of the 
ABAB gauge we typically used a cutoff of A = 10, where 
\Qx\> \Qy\ < A and Qi are the coefficients of the recipro- 
cal lattice primitive vectors. For the AABB gauge the 
unit cell is non-square, so we took 2|Q X |, \Q y \ < A, typ- 
ically with A = 8. It was not feasible to push A much 
greater than 20. Also, for any fixed A, we expect the 
numerical calculations to become less reliable as ajj is 
increased, since raising the anisotropy introduces more 
and more low-energy states that do not fall within the 
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cutoff and are improperly thrown out. Although the ac- 
cessible range of A and old is not sufficient to provide 
detailed information on the true asymptotic behavior as 
old —* oo or the convergence of the numerics, we believe 
it should be sufficient to extract certain important qual- 
itative features. 




FIG. 4. Plot of ln(v-(l)/v-(ao)) versus the anisotropy 
old, comparing data from the 2d numerics with the Id model. 
The results shown are in the AABB gauge, with tilt angles 
giving distinct columns in the Id model. The hollow shapes 
are numerical results for 6 = and 8 — arctan 2. The lines 
are the asymptotic results from the Id model and have slope 

7T 2 /(16dS). 

Our results for V-(old) at fixed 8 show good qualita- 
tive agreement with the predictions of the Id model. For 
tilt angles and gauges with distinct A and B columns 
we see roughly exponential decrease of t>_ similar to that 
calculated analytically (Fig. 4). We believe we do not 
see quantitative agreement because we cannot access the 
asymptotic regime numerically. When there are over- 
lapping columns we see similar agreement with the ana- 
lytical results; t>_ is roughly proportional to l/a D with 
the constant of proportionality closer to unity when the 
columns are farther apart (Fig. 5). This is consistent 
with our understanding of the Id model; at greater col- 
umn separation d v is smaller, and the quasiparticle wave- 
functions do not have to be smeared out as much for them 
to be completely delocalized in the y-direction. Numer- 
ical calculations of v + show only small variations, typi- 
cally by at most a factor of order unity as old ranges from 
1 to 8. Therefore the renormalized anisotropy v + /v- be- 
haves similarly to 1 /v— as a function of old ■ 

We have seen that, for large ajj, tilting of the vortex 
and crystal lattices in the Franz- Tesanovic framework can 
be understood in terms of a one-dimensional model. For 
tilt angles and singular gauges with distinct columns of 
A and B vortices, w_ is exponentially small in etc, while 
when the columns overlap it decreases algebraically. In 
both these cases v + does not vary significantly with old ■ 



These results depend dramatically on the choice of singu- 
lar gauge, and provide a further indication that the FT 
transformation must somehow be regularized if it is to 
provide a consistent account of the physics. While pro- 
posals for the kind of regularization necessary have been 
discussed in the literature 13 ' 15 ' 16 , a detailed resolution 
of these issues remains an important problem for future 
work. 




FIG. 5. Plot of 1/v- versus an comparing 2d numerical 
data in the AABB gauge with the Id model. The tilt angles 
are chosen to give overlapping columns of vortices in the Id 
model. The shapes connected by dotted lines are 2d numeri- 
cal data for different tilt angles, and the solid line is simply a 
plot of 1/v- — an, the exact Id result. 
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